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Abstract. The nucleon structure and the nuclear force are investigated in the context of the non- 
perturbative path-integral method of hadronization. Starting from a microscopic quark-diquark 
model, the nucleon is generated as a relativistic bound state and an effective chiral meson-nucleon 
Lagrangian is derived. Many of the nucleon physical properties are studied using a theory of at most 
two free parameters. 


INTRODUCTION 

The central challenge in nuclear physics remains to understand the origin and nature 
of the nuclear force due to our inability to solve quantum chromodynamics (QCD), the 
fundamental theory for the strong interactions. The basic problem of QCD is that its 
fundamental degrees of freedom, quarks and gluons, are not the observable baryon and 
meson states. Thus bridging the missing link between the fundamental and observable 
degrees of freedom stands as one of the stark challenges of nuclear/elementary particle 
physics today. Although we do have an ab initio approach to solve this problem, that is 
lattice QCD, this endeavor is still miles away from achieving such a goal. For the time 
being, we have no alternative but to resort to effective non-perturbative approaches of 
which this study is one. 

This presentation describes our work [1] in addressing this missing link by deriving 
a chiral meson-nucleon Lagrangian from a microscopic model of quarks and diquarks 
using the path-integral method of hadronization. Chiral symmetry and its spontaneous 
breaking have proven to be key concepts in understanding meson and baryon structure 
and many features of the nuclear force. Therefore, we start from a QCD-based chiral 
effective field theory, the Nambu Jona-Lasinio (NJL) model that accommodates most 
of QCD symmetries [2]. Then, the nucleon is described as quark-diquark correlations. 
This assertion is vindicated by a mounting experimental evidence that diquarks play a 
dynamical role in hadrons [3]. Using the path integral hadronization, we calculate nu¬ 
cleon properties and derive an effective chiral meson-nucleon Lagrangian of the quan¬ 
tum hadrodynamics (QHD) type [4] that describes the rich meson-nucleon interactions 
in a fully covariant and chirally symmetric formalism. 

While this program is applied to the case of nucleons and mesons, it is certainly 
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of general nature and can possibly be applied prolifically to yield other baryons and 
their interactions. Moreover, the idea of using path-integral techniques to transform a 
Lagrangian from its fundamental to its composite degrees of freedom is a powerful 
concept in physics of immense impact and utility. As a matter of fact, the authors 
of Ref. [5] have recently invoked such techniques in their study of high-temperature 
superconductivity. They succeeded in doing so by converting a model of strongly- 
correlated electrons into an effective U(l) gauge field theory in terms of composite fields. 

The use of hadronization has been introduced in Ref. [6, 7]. Consequently, the authors 
of Ref. [8] attempted to construct an effective Lagrangian for the nucleon using only 
scalar diquarks. In the present work, we extend their work by deriving the structure 
using both axial-vector and scalar diquarks and we employ, as opposed to Ref. [8], a 
gauge-invariant regularization scheme throughout our analysis. Furthermore, we verify 
the Ward identity and the Goldberger-Treiman relation and present a full numerical study 
of various nucleon observables for the case of scalar diquarks drawing special attention 
to the role of an intrinsic diquark form factor. Hence, this work is the first calculation 
of an extensive set of nucleon observables using path-integral hadronization since the 
introduction of the idea more than ten years ago. 


FORMALISM 

We start from an NJL Lagrangian satisfying SU(2) L x SU(2)^ chiral symmetry: 

%njl = <?(*£ - m o)q + 2 [(w) 2 + (^‘y 5 ^) 2 ] , (i) 

where q is the current quark field, f are the Pauli matrices, G is the NJL coupling 
constant, and m 0 is the current quark mass which explicitly breaks chiral symmetry. 
The color and flavor indices are suppressed for brevity. 

By introducing composite scalar (o ~ qq) and pseudoscalar (k ~ qiy 5 Tq) fields 
through the Hubbard-Stratonovich transformation [9], we can rewrite the NJL lagrangian 
into a semi-bosonized Lagrangian ^f^ JL — [q(i9 — a — iy s T ■ 7i)q — jq(o 2 + K 2 )], 
where we have absorbed the bare quark mass m Q into the sigma field a. Further, we 
transform the meson fields according to the non-linear parameterization [<r,7r] —> [d / ,T>]: 

a + iy 5 z ■ K = (m q + o') . Here, F n = 93 MeV is the pion decay constant and 

m q = (<x) 0 is the constituent quark mass which is fixed through a gap equation in the 
meson sector [2]. 

As a consequence of considering Lorentz structure, there are five types of possible 
diquark qq correlations. These are the scalar, pseudo-scalar, vector, axial-vector, and 
tensor diquarks. For each of these Lorentz qq formations, we have also an isoscalar 
and isovector diquarks. Using permutation symmetry and Fierz transformation, we have 
verified an earlier assertion [10] that only two diquark formations are independent for 
the nucleon if its field is to be written as a local operator of three quarks. Hence, we 
introduce />' as an axial-vector isovector diquark field and D as a scalar isoscalar one. 

Next we introduce a quark-diquark interaction term in such a manner to generate the 
nucleon as a linear combination of axial-vector and scalar diquarks. It is convenient 
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quark field defined by % = e F n - q. By introducing electromagnetic interactions and 
batching the semi-bosonized NJL Lagrangian, the diquark contributions and the quark- 
diquark interaction term, we obtain the following Lagrangian as our microscopic model: 
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Here 0 is a mixing angle for the two diquark contributions, G is the quark-diquark 
coupling constant and Q q , Q s and Q A are the quark and diquark charges while S Q , 
A 0 and A 0 ((y are the free quark and diquark propagators (Notice that &{Q 2 SA ) terms 
are discarded in Eq. (3)). The matrix contains all interaction vertices of the quark 
field with meson and electroweak fields (weak part is not shown), where the vector 

and the axial vector srffi fields are defined through the Cartan decomposition (c, = y-)\ 

exp^-iy 5 ?-!) ^ exp^T-f) = 

Subsequently, we introduce collective nucleon fields (B ~ sin0D v ■ fy v y 5 ^ + 
cos 9D%) through another Hubbard-Stratonovich transformation. At this point only 
quarks and diquarks are dynamical fields with kinetic terms while the meson and 
nucleon fields are merely auxiliary ones. By integrating over the quark and then over the 
diquark fields, we obtain a meson-nucleon effective Lagrangian. Accordingly, we have 
“hadronized” the microscopic theory by producing the dynamical meson (bosonization) 
and nucleon (fermionization) fields. 

Thereupon, we arrive at a compact Lagrangian given by 

^f eff = SJf sb — ■^—(o , + in q ) 2 — i tr lnS 1 1 -L BB + ztrln(l — □) + 

2 G G 

i trln(l — A 0 EMInt) + ztrln(l — A 0 EM Int) . 
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FIGURE 1. The self-energy diagram which generates the nucleon kinetic and mass terms and produces 
the mass equation that determines the nucleon mass. 


Here the trace is over color, flavor and Lorentz indices while the “EM Int” label stands 
for the diquark electromagnetic interaction terms. Furthermore, 
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This effective Lagrangian contains plenty of rich physics: kinetic and mass terms for nu¬ 
cleons and mesons together with a multitude of interaction terms of mesons, nucleons, 
and electroweak gauge bosons. Nonetheless, the most desired part of the Lagrangian is 
the prized chiral meson-nucleon interaction and nucleon-nucleon vertices which delin¬ 
eates the nuclear force. 

In order to explore the physics of the nucleon sector, we take the lead¬ 
ing term in the loop and derivative expansion of ztrln(l — □) — T BB —> - 

J d 4 xd 4 yB(x) E(jc,y) + ^8(x~y) B(y), which is nothing but the nucleon self¬ 
energy (see Fig. 1). The Fourier transform of E is then decomposed as E(p) = 
Up 2 ) + i } ^v(p 2 )- This leading term generates dynamically the nucleon mass M B 
which is extracted as the pole of the propagator: 


l + U M l) + M r UM 2 b) = 0 . 
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Thus, near the mass shell the inverse nucleon propagator takes the form: 
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where Z is the wave-function renormalization constant (B = \fz 0 ren ). Evidently, the 
nucleon has finally acquired the desired status as a dynamical degree of freedom in the 
problem. 

In computing the various Feynman diagrams in the problem, we encounter divergent 
integrals that must be regularized. Several regularization schemes were attempted. We 
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FIGURE 2. The Feynman diagrams for the electromagnetic coupling which generate the electromag¬ 
netic vertex of the nucleon and subsequently are used to test the validity of the Ward identity. 


found that the most suitable scheme is the PV technique which we have adopted as the 
standard method in this work. Accordingly, we verified the Ward-Takahashi identity by 
computing the electromagnetic vertex shown in Fig. 2. As a matter of principle, the PV 
mass in the nucleon sector can be different from the NJL cut-off arising in the meson 
sector [2]. Nonetheless, to minimize the number of free parameters, we elected to equate 
them. It is noteworthy here that all observables were found to be very insensitive to the 
value of the PV mass upholding the futility of using it as a free parameter. 

We have also verified one of the chiral symmetry relations, the Goldberger-Treiman 
relation, by computing both the axial-vector coupling constant g A and the pion nucleon 
coupling constant g nNN - In the hadronization formalism, this relation emerges naturally 
intact as opposed to large violations in the Bethe-Salpeter equation approach [11] 


NUMERICAL STUDY 

Having derived the structure of the Lagrangian, we proceed to generate numerical results 
using only scalar diquarks (6 = 0 in Eq. (2) and (5)), thereby admitting the possibility 
of an intrinsic diquark form factor (IDFF). Including only scalar diquarks is not out of 
place as many recent studies using such diquarks have reported good results for most 
of the nucleon observables [12, 13, 11]. Moreover, there are strong indications of large 
scalar diquark dominance in the nucleon [14]. 

The parameters G and the cut-off A are fixed to yield the constituent quark mass 
and the pion decay constant [2]. The diquark masses are also determined in the NJF 
model [15]. This leaves us with only one free parameter in our model: the quark- 
diquark coupling constant G which is fixed to determine the nucleon mass of 0.94 GeV 
through the mass equation (6). Thus the basic quantities in our model are the constituent 
quark mass m q — .390 GeV, the scalar diquark mass M D = .600 GeV, the quark-diquark 
coupling constant G = 159.1 GeV -1 , and the Pauli-Villars mass A = .600 GeV. We 
obtain a binding energy of AE hjn = m q + M D — M B = 50 MeV, suggesting a loosely 
bound state for the nucleon. 

Tab. 1 displays our predictions for some of the static properties of the nucleon. For 
the nucleon magnetic moments, our treatment predicts a number that is two-third of 
the experimental value for the proton and one-half of that for the neutron. This is 
not a surprising result as we have not included the axial-vector diquark in the present 
calculation. The predicted value for the axial-vector coupling g A of 0.87 is less than the 
experimental one of 1.26 indicating here also the importance of the axial-vector diquark. 


TABLE 1 . Some of the nucleon static properties as predicted in the present calculation using 
the intrinsic diquark form factor (IDFF) or without it. Experimental values are taken from 
Ref. [16, 17], 
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FIGURE 3. The nucleon electric form factors with and without the intrinsic diquark form factor (IDFF) 
along with its quark and diquark contribution. The left panel shows the proton results while the right panel 
displays those for the neutron. Experimental data can be found in Ref. [19, 20]. 


The nucleon size is nicely well-produced in our model as the electric and magnetic 
radii for the nucleon are very close to the experimental measurements. The negative 
charge radius of the neutron has been suggested as an indication of a scalar diquark 
clustering in the nucleon [18] and our treatment dynamically manifests this assertion. 
These numbers point to a physical picture of a “heavy” diquark at the center with a 
quark rotating around it. The extended size of the diquark contributes a positive value 
of about 0.10 fm 2 for the nucleon electric radii. As expected, the IDFF has virtually no 
effect on the magnetic radii as the scalar diquark has a negligible contribution to the 
magnetic form factors. 
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FIGURE 4. The nucleon magnetic form factors with and without the intrinsic diquark form factor 
(IDFF) along with its quark and diquark contribution. The left panel shows the proton results while the 
right panel displays those for the neutron. Experimental data can be found in Ref. [19, 20]. 


Next we calculate the nucleon form factors. The left panel of Fig. 3 displays the proton 
form factor with and without the IDFF along with its quark and diquark contribution. Our 
treatment produces beautifully this observable. It is evident here that the IDFF [21] plays 
an important role specially at large values of momentum transfer ( Q 2 = —q 2 where r/' 
is the momentum transfer). The neutron electric form factor tells a similar story (right 
panel). Clearly, the quark contribution is negative in value (d-quark) and thus cancels 
much of the diquark contribution leading to a small form factor. It is noteworthy here 
that the neutron form factor is a potent test of any treatment as it is a delicate cancellation 
of two large contributions [13]. Saliently, the cancellation is naturally produced in our 
study. 

In Fig. 4 we present the nucleon magnetic form factor as calculated with or without 
the intrinsic diquark form factor. Unmistakably, the scalar diquark contribution is vir¬ 
tually vanishing due to the lack of an intrinsic spin. Nevertheless, there is a very small 
contribution due to a small orbital angular-momentum effect in the bound quark-diquark 
system. A comparison with experimental data suggests the need for the axial-vector di¬ 
quark, which does have an intrinsic spin, to supplement the quark contributions and to 
provide the missing strengths for the magnetic form factors. 










CONCLUSIONS 


In conclusion, we have tackled the nucleon structure and the challenging problem of 
understanding the origin and nature of the nuclear force by deriving a meson-nucleon 
Lagrangian using the path-integral method of hadronization. The treatment produced a 
remarkable agreement with experimental data for the nucleon size and its electric form 
factors, while our calculations show missing strengths for the magnetic properties and 
g A . The discrepancy is likely due to the absence of the axial-vector diquark in the present 
numerical study. This presentation describes the first work in our program of using path- 
integral hadronization to study baryon structure and the nuclear force. Deriving this force 
provides nuclear physics with a solution to the stigma of no fundamental foundation that 
has tarnished its image for decades. 
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